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Abstract 



(-H We investigate the role of centrifugal acceleration of electrons in producing the very high energy (VHE) radiation from the BL Lac 
O pbiect lES 0806+524, recently detected by VERITAS. The efficiency of the inverse Compton scattering (ICS) of the accretion disk 
I thermal photons against rotationally accelerated electrons is examined. By studying the dynamics of centrifugahy induced outflows 
2 a-nd by taking into account a cooling process due to the ICS, we estimate the maximum attainable Lorentz factors of particles 
and derive corresponding energetic characteristics of the emission. Examining physically reasonable parameters, by considering the 
^ Iiarrow interval of inclination angles (0.7° — 0.95°) of magnetic field lines with respect to the rotation axis, it is shown that the 
I icentrifugally accelerated electrons may lead to the observational pattern of the VHE emission, if the density of electrons is in a 
certain interval. 
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1. Introduction 

■ In physics of active galactic nuclei (AGNs) one of the 
major problems is related to the understanding of origin of 
the high energy radiation. One prominent class of AGNs 
is the so-called BL Lac objects - supermassive black holes 
characterized by rapid and large amplitude flux variability. 
By using the radio observations from the Green Bank 91-m 
telescope [I], the AGN, lES 0806+524, was identified as a 
BL Lac object [2]. 

Recently, by VERITAS was found that the blazar, lES 
0806+524 reveals VHE spectra in the TeV domain [3 . Ac- 
cording to the standard model of BL Lacs, VHE radia- 
tion originates from the ISC of soft photons against ultra- 
relativistic electrons [4 5 . However, the origin of efhcient 
acceleration of particles up to highly relativistic energies 
still remains uncertain and needs to be revealed. Proposed 
mechanisms based on the Fermi-type acceleration process 
[B] may be applied successfully for the TeV emission, only, 
if the initial Lorentz factors of electrons are considerably 
high (7 > 102) [7. 

It is clear that in the rotating magnetospheres (the inner- 
most region of AGN jets and pulsar magnetospheres) the 
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centrifugal effect should play a significant role in the overall 
dynamics of corresponding plasmas. For example, the ro- 
tationally driven parametric plasma instabilities have been 
studied for pulsars 8 9 and AGNs |10|11) respectively, and 
was shown, that under certain conditions, the relativistic 
effects of rotation may efficiently induce plasma instabili- 
ties, parametrically pumping the rotational energy directly 
into the plasma waves. The centrifugahy induced outflows 
have been discussed in a series of works. Blandford & Payne 
in the pioneering paper [12j considered the angular mo- 
mentum and energy pumping process from the accretion 
disk, emphasizing a special role of the centrifugal force in 
dynamical processes governing the acceleration of plasmas. 
It was shown that the outflows from accretion disks oc- 
curred if the magnetic held lines are inclined at a certain 
angle to the equatorial plane of the disk. In the context of 
studying the nonthermal radiation from pulsars, the cen- 
trifugal effect has been examined in [13114115] , where the 
curvature emission of accelerated particles was studied. By 
applying the similar approach, Gangadhara & Lesch con- 
sidered the role of centrifugal acceleration on the energetics 
of electrons moving along the magnetic held lines of spin- 
ning AGNs [16j . This work was reconsidered in a series of 
papers [7|17|18j and the method was applied to a special 
class of AGNs - TeV AGNs. It was shown that considera- 



tion of straight field lines is a good approximation and was 
found that the centrifugal force may accelerate electrons 
up to very high Lorentz factors (~ 10^) providing the TeV 
energy emission via the ICS. 

In the present paper we investigate a role of rotational 
effects in the VHE flare from the blazar lES 0806+524, by 
applying the method of centrifugal outflows, developed in 
[7117118119] . We show that, for a certain set of parameters, 
due to the ICS in the Thomson regime, photons, when up- 
scattered against centrifugally accelerated ultra-relativistic 
electrons, produce the VHE radiation in the TeV domain. 
We show that a resulting luminosity output is in a good 
agreement with the observed data. 

The paper is arranged as follows. In Sj2]we consider our 
model and derive expressions of the luminosity output and 
the energy of photons respectively. In Sj3] we present the 
results for the blazar lES 0806+524 and in ^we summarize 
our results. 

2. Main consideration 

Let us consider the typical parameters of lES 0806+524: 
the black hole mass, Mbh ~ 5 x IO^Mq [H], {Mq is the so- 
lar mass) and the bolometric luminosity, L w 7 x lO'^'^erg/s 
[22] . We examine particles originating from the accretion 
disc at the distance ^ 10 x Rg from the central object, 
where Rg = 2GMbh / is the gravitational radius of the 
black hole. Then, by taking the value of the equipartition 
magnetic field. 



(1) 



into account, one can show that for typical parameters, r « 
10xi?g,n G (0.0001 — l)cTO~^, 7o ~ 1, the value of the ratio, 
/ ^Qmnc^ , is in the following interval ^ 10^ — 10^'^ (70, n 
and m are electrons' initial Lorentz factor, the density and 
the rest mass respectively). Therefore, the magnetic field 
energy density exceeds the plasma energy density by many 
orders of magnitude, which indicates that the plasma co- 
rotates with the angular velocity. 
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corresponding to the Keplerian motion at tq ~ 10 x _Rg. 

We see that due to the frozen-in condition the parti- 
cles follow the co-rotating magnetic field lines and accel- 
erate centrifugally. Therefore, it is reasonable to consider 
dynamics of the electron, sliding along the rotating mag- 
netic field lines. We apply the method developed for AGNs 
in [17118] and assume that the straight field lines co-rotate. 
Then, if we take an angle a between the magnetic field, 
B, and the angular velocity of rotation, uj, into account, 
after the transformation of coordinates: x = rsinacosujt, 
y = rsinasinojt and z = rcosa of the Minkowskian metric, 
ds"^ = rjfj.udx'^dx^ (77^1/ = diag{—l, and x'^ = 



[ct] X, y, z)), the metric in the co-moving frame of reference 
is given by [17119] 



For the equation of motion we get: 

d dC _ dC_ 
dx Q ^dsi^'^ dxt^' 

1 _ dxf" dx" 



gf,^ = diag <- [1 
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where 



~ a; sin a, 
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Then, by taking the four velocity identity, gafs = 
— 1, into account, one can derive from Eq. ([4]) the radial 
equation of motion fT9\ : 



dt^ 
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Solving Eq. ([7|), it is straightforward to show that the 
Lorentz factor of the particle changes radially as ^ : 
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where 
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To and Wo are the initial position and the initial radial ve- 
locity of the particle, respectively and Ric is the radius of 
the light cylinder - a hypothetical zone, where the linear 
velocity of rigid rotation exactly equals the speed of light, c. 

As is clear from Eq. ([5]), in due course of time the Lorentz 
factors of electrons become very high in the vicinity of the 
light cylinder (r ~ Ric)- On the other hand, it is clear that 
acceleration lasts until the electron encounters a photon, 
which in turn inevitably limits the Lorentz factor of the par- 
ticle. During the ICS an electron will lose energy, whereas 
a photon will gain energy. This mechanism is characterized 
by the so-called cooling timescale [2D] 



cool 



= 3 X 
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where Urad — L/Aircr'^ is the energy density of the radia- 
tion. The acceleration process is characterized by the accel- 
eration timescale, tacc = ^/(d-y/dt), which after applying 
Eq. ([5]) can be presented by 



tacc — 
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Generally speaking, initially the electrons accelerate, but 
in due course of time the role of the inverse Compton losses 
increase and the acceleration becomes less efficient. The 
maximum energy attainable by electrons is achieved at a 
moment when the energy gain is balanced by the energy 
losses due to ICS. Mathematically this means that the fol- 
lowing condition tacc ~ tcooi has to be satisfied. After ap- 
plying Eqs. (|9ll0p the aforementioned condition leads to 
the expression of the maximum Lorcntz factor j2] 
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where Ri ~ Rid sin a. If electrons with such high kinetic 
energies encounter soft photons having energy, tg, then, 
photons' energy after scattering is given by 



e = 7„ 



(12) 



As we have already mentioned, the particles reach maxi- 
mum kinetic energy almost on the LC surface. Let us as- 
sume that a layer where the ICS takes place and the high 
energy photons are produced has a thickness, Ar. Then, 
for the corresponding infinitesimal volume of a cylindrical 
layer we get: 
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If we take a single particle Thomson power 

Pic ^ (TTC^^Urad , 
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into account, then the total power emitted from the radia- 
tion zone can be expressed as follows 
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da, (15) 



where ax ~ 6.65 x 10~ cm is the Thomson cross-section. 

3. Discussion 

According to the observations of VERITAS, per- 
formed from November 2006 to April 2008, the blazar 
1ES0806-I-524, has the VHE emission with the intrinsic 
photon spectrum, Ti = (2.8 ± 0.5), and the intrinsic inte- 
gral flux, J^>0 3^,^, = (4.4 ±1.1) X 10-i2cm-2s-i ^3J. Then 
the intrinsic differential photon spectrum can be easily 
restored by the integral flux 
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where Fo « (11.8 ±2.9) x IQ-^'^ cm-'^s-^. 



The corresponding luminosity above O.STeV up to ITeV 
can be estimated as follows: 



AS l^ede, 
de 



(17) 



where ei — Q.STeV and 62 = ITeV. We assume that 
the high energy emission originates from the jet, hav- 
ing an opening angle, 2q!,„ (where a^n > 0^2 )• AS* « 
ttD^ (sin^ a2 — sin^ ai) and D sa 630Mpc is the distance 
to the blazar. 

By taking the parameters into account, one can see from 
Eq. pT|) that the luminosity in the energy interval (0.3 — 
l)TeV is given by 



LvHE ~ (3-7 ± 1.3) X lO''^ (sin^ 0^2 - sin^ ai) 



erg 



(18) 



According to the standard theory, it is well known that 
the accretion disks thermally radiate and the corresponding 
temperature is expressed as in the following way 124] : 



T 



3.1 X 10«Q2/5^^-l/5rf-3/5 (^^ _ ^-1/2^ ^ ^ (^9) 



where 



Q = 



M yr 
3M8 yr- 



(20) 



is the dimensionless mass accretion rate, Ms = Ms^f/IO^M© 
and d = r^/3Rg. 

For the given luminosity, the mass accretion rate can be 
estimated as: 
L 



M w , (21) 

then, combining Eqs. (|19l2ip one can show that energy, 
eg — kT, of accretion disc's thermal photons emitted in the 
area from r* = 15 x Rg to r* = Ric is of order ~ lOeT^. 
Therefore, as we see from Eq. ()12|) . for producing energies 
from thousands of GeV to TeV domain, one requires very 
high Lorentz factors (1 — 3) x 10^. 

One can see that the aforementioned values of Lorentz 
factors are achieved for very low inclination angles. In Fig.[l] 
we show jmax as a function of the inclination angle. The set 
of parameters is vq = 0.4c, ro = 10 x Rg, u = 4.5 x 10^®s~^ 
and L = 7 X lO'^'^erg/s. As is clear from the figure, the 
electrons reach high values of the Lorentz factor for small 
angles, 0.7° — 0.95°. This is a natural result because, one 
can straightforwardly show from Eq. (jlip . that 'jmaxio:) 
behaves as 1/ sin^ a and therefore, provides higher kinetic 
energies for lower inclinations. 

The present model is based on an assumption that max- 
imum kinetic energy of particles is determined by the bal- 
^ _ ance of energy gain due to the acceleration and energy 
losses due to the ICS. Generally speaking, this approach is 
valid only if the energy losses is dominated by the ICS. On 
(16) the other hand, apart from the inverse Compton scatter- 
ing, also the curvature radiation could impose significant 
limitations [25l . 
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Fig. 1. The behaviour of the maximum Lorentz factors of electrons 
versus the inchnation angle. The set of parameters is vq = 0.4c, 



Fig. 2. The dependence of the energy of photons on the inclina- 



ro = 10 X Rg 



■■ 4.5 X IQ-^s-l and L = 7 X 



The centrifugal acceleration mainly happens close to the 
light cylinder, and since the power of a single particle cur- 
vature radiation behaves as ^ 'y^, one has to check the 
constraint imposed by this mechanism on relativistic par- 
ticle dynamics. A total power radiated by a single particle 
is given by 



P = 2 e^^4 



(22) 



where by Rc we denote the curvature radius. Then, the 
timescale of curvature emission can be defined by the fol- 
lowing way: 



^rac 



(23) 



To find the limitation imposed on the maximum Lorentz 
factor let us note that electrons initially accelerate effi- 
ciently, and this process lasts until the energy gain is bal- 
anced by the curvature losses. This happens when tacc ~ 
tc- By taking Eqs. (|10l23p into account and assuming Rc ^ 
Ric, it is straightforward to show 
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From Eqs. (|ll|24p we see that for 70 ~ 1 one has the 
following inequality j^^x ^ Jmax- This indicates that the 
curvature radiation does not impose a significant limita- 
tion on the maximum attainable Lorentz factors. Therefore 
we conclude that maximum attainable kinetic energies are 
determined only by the ICS. 

As is clear from Eq. the VHE emission (0.3reT^ - 
ITeV) detected by VERITAS can be achieved for the fol- 
lowing interval of Lorentz factors (1.5 — 2.8) x 10''^. Equa- 
tion (|lip shows that when the inclination angles are in the 
following range a G (0.7° — 0.95°), electrons attain Lorentz 



tion angle. The set of parameters is vq = 0.4c, ro 
u) = 4.5x 10~6<(-i and L = 7 X 10**erg/s. 



10 X R„ 



factors from the aforementioned region of values (see Fig. 
[I]). Wo = 0.6c, ro = 10 X Rg, = 4.5 x lO^^s^^ and L = 
7 X 

In Fig. [2] we show the behavior of the emission energy 
versus the inclination angle. As is clear from the figure, ac- 
celeration of electrons inside the region, 0.7° < a" < 0.95°, 
provides the photon energies from 0.3Tev up to ITeV. 

For finding the luminosity output above 0.3rel^, we have 
to integrate Eq. from ai sa 0.7° to 02 ~ 0.95°. Let 
us consider the following set of parameters: n = 1.1 x 
10~^cm~^, Vq — 0.4c, ro = 10 X Rg and r* = 15 x Rg, 
then, the integration in Eq. gets the VHE luminosity 
output of order 4.8 x lO^^erg/s. From Eq. we see that 
for ai < a < a2, the observed VHE luminosity is in the 
range ~ (4.6 ±1.6) x lO^^erg/s, therefore, our results are 
in a good agreement with the observed data. 

For plotting our graphs and getting the results we used 
the parameters with the best fitting, although the values 
from the following ranges are also applicable: ro/Rg E [10 — 
20], Vq e [0- 0.7c], n e [0.7 - 1.4] x IQ-^cm-^^ r*/Rg e 
[10-20]. 

On the other hand, the high energy photons may un- 
dergo the 77 absorbtion. It is well known that gamma rays 
interact most effectively with the background photons of 
energy [26j 



^(mc2)2 TeV 
eb = 4^ — « 1 eV 
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and the corresponding cross section has a peak at ctq R 
<7t/^- The optical depth of high energy photons, then be 
comes 

Ric 



T = — = niROoRi 



where A is the mean-free path of infrared photons. 



niR = 
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(26) 



(27) 



A 



is the corresponding photon density and L{eh) - the infrared 
himinosity. After substituting Eq. (l27l) into Eq. ((26)) . one 
can derive an expression of the optical depth of high energy 
photons [21] 

where L^dd ~ 6.5 x lO^^erg/s is the Eddington luminosity 
of lES 0806+524. 

As we see from the figures, for producing radiation in the 
ITeV domain, one has to accelerate the electrons up to 7 « 
2.8 X 10^. For lES 0806+524 no infrared data are published 
so far, on the other hand, since the TeV emission is detected, 
therefore the infrared luminosity of lES 0806+524 must be 
less than 5.6 x lO-^L^dd « 3.6 x 10^°erg/s (see Eq. ([28]) ). 
One has to note that the centrifugal acceleration leads to 
the TeV variability timescale of the order ~ {l~2)days [18 . 
It is worth noting that TeV blazars exhibit the variability 
on hour to minute timescales, but this particular feature is 
not detected for lES 0806+524. 

As our investigation shows, the centrifugally accelerated 
outflows may provide the detected VHE emission of lES 
0806+524 via the ICS if the parameters are chosen appro- 
priately. From the aforementioned set of parameters very 
important one is the density of relativistic electrons, which 
has to be in the following interval [0.7 — 1.4] x 10~^cm~'^ 
in order for the centrifugal acceleration to explain the de- 
tected high energy emission. This we consider as a certain 
test to check if the mentioned mechanism is feasible. 

4. Summary 

(i) For explaining the observed TeV energy radiation 
from lES 0806+624 detected by VERITAS, we have 
considered the inverse Compton scattering of disk 
thermal photons against centrifugally accelerated 
ultra-high energy electrons. 

(ii) We have shown that due to very strong magnetic field, 
the electrons are in the frozen-in condition, which 
leads to the co- rotation of particles. Due to the co- 
rotation, electrons centrifugally accelerate almost up 
to the light cylinder surface, and in the nearby zone 
of it the electrons upscatter against soft thermal pho- 
tons, causing the limitation of particles' kinetic en- 
ergy. We also have shown that the role of the curva- 
ture radiation in limiting the maximum kinetic en- 
ergy is negligible with respect to the inverse Compton 
losses. 

(iii) We considered the 77 absorbtion of high energy pho- 
tons in the background field of soft infrared photons. 
From the observationally evident fact that the TeV 
radiation escapes the central source, we estimated the 
maximum value of infrared luminosity which has not 
been detected so far. 

(iv) We have found that for physically reasonable param- 
eters the ICS occurs in the Thomson regime. It has 



been shown that for the following interval of inclina- 
tion angles, 0.7° — 0.95°, with the best fitting parame- 
ters the resulting emission energies, {O.STeV— ITeV) 
and the luminosity output, lO'^^erg/s are in a good 
agreement with the observed data. On the other hand, 
if the parameters are chosen in certain physically rea- 
sonable intervals, the high energy emission is also pos- 
sible. Therefore we offer the following test: if one in- 
directly measures the density of the relativistic elec- 
trons, and finds its value to be in the following range, 
[0.7 — 1.4] X lO^'^cTO^'^, then the centrifugal acceler- 
ation is a feasible mechanism in producing the TeV 
photons via the ICS. 
In the paper we have done several approximations. The 
first limitation is that we studied the straight magnetic field 
lines, although, especially in the very vicinity of the light 
cylinder the curvature of field lines becomes significant. 
Therefore, the generalization of the present approach will 
be the next objective of our future work. 

The next approximation concerns the fact that accord- 
ing to our model magnetic field is not infiuenced by plasma 
kinematics. On the other hand, for real astrophysical sce- 
narios, it is obvious that plasmas may undergo the overall 
configuration of the magnetic field. For this reason, it is 
very important to generalize the approach presented here 
and see how the collective phenomena change the results. 
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